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I. INTRODUCTION 

The strong force — one of the four fundamental interac- 
tions of nature along with gravity, electromagnetism and the 
weak force — is the force that holds together protons and 
neutrons in the nucleus. The strong interaction is described 
by Quantum Chromodynamics (QCD) |1]. This description 
is based on a model of elementary particles — the quarks 
— possessing "color charge" and interacting through the ex- 
change of gauge fields — the gluons (equivalent to the pho- 
tons in the electromagnetic interaction). QCD is a quantum 
field theory, with local SU(3) gauge symmetry, correspond- 
ing to three possible colors. The fact that the gauge group 
of QCD is non-Abelian implies that the gluons possess color 
charge and therefore interact with each other, as opposed to 
the photons. The only parameters of the theory are the masses 
of the various types (called "flavors") of quarks considered 
and the value of the strong coupling constant. 

A unique feature of the strong force is that the particles 
that feel it directly — quarks and gluons — are completely 
hidden from us, i.e. they are never observed as free particles. 
This phenomenon, known as confinement, makes QCD much 
harder to handle theoretically than the theories describing the 
weak and electromagnetic forces. Indeed, the coupling con- 
stant oc.s of the strong interaction becomes negligible only in 
the limit of small distances, or equivalently in the limit of high 
energy or momentum. This property is called asymptotic free- 
dom. At larger distances (i.e. smaller energies) there is an in- 
crease in the intensity of the interaction and it is believed that 
the force of attraction between quarks is constant, for suffi- 
ciently large distances, determining the confinement of quarks 
and gluons inside the hadrons. The fact that a. v is not negli- 
gible at low energies makes the study of important phenom- 
ena such as the mechanism of quark and gluon confinement, 
the hadron mass spectrum and the deconfining transition at 
finite temperature inaccessible to calculations using perturba- 
tion theory. These phenomena must therefore be studied in a 
nonperturbative way. 

The nonperturbative study of QCD is possible in the lattice 
formulation of the theory, introduced by Wilson in 1974 Q], 
This formulation offers a convenient nonperturbative regular- 
ization, preserving the theory's gauge invariance. The essen- 
tial ingredients for the lattice formulation are: 1) path-integral 
quantization, 2) continuation to imaginary or Euclidean time 



and 3) lattice regularization (given by the discretization of 
space-time). The combination of the first two ingredients 
makes the theory equivalent to a model in classical statistical 
mechanics: indeed, in Euclidean space a path integral for the 
quantum theory is equivalent to a thermodynamic average for 
the corresponding statistical mechanical system. For QCD, 
the square of the bare coupling constant go of the field theory 
corresponds directly to the temperature 1/p of the statistical 
mechanical model. 

The third ingredient — the lattice discretization — repre- 
sents an ultraviolet regularization. In fact, the lattice spac- 
ing a corresponds to a high-momentum cutoff, since momenta 
higher than ~ 1/a cannot be represented on the lattice. In this 
way the ultraviolet divergences, appearing in the calculation 
of physical quantities, are suppressed and the theory is well 
defined. Of course, in order to recover the continuum-space 
theory we must take the limit a — ► 0. In this process it is nec- 
essary to "tune" the bare parameters of the theory in such a 
way that physical quantities converge to finite values, which 
can then be compared to experiment. In particular, in the limit 
a — > 0, a correlation length Z, measured in units of the lattice 
spacing, i.e. %/a, must go to infinity. In other words, the lat- 
tice theory considered must approach a critical point, i.e. a 
second-order phase transition. Thus, the study of the contin- 
uum limit in quantum field theories on the lattice is analo- 
gous to the study of critical phenomena in statistical mechan- 
ics. The correspondence between Euclidean field theories and 
classical statistical mechanics allows the application of usual 
statistical-physics methods to the study of QCD. In particu- 
lar, one may perform numerical simulations by Monte Carlo 
methods, which are based on a stochastic description of the 
systems considered |3]. 

Despite the similarity of the methods, the Monte Carlo sim- 
ulation of gauge theories is much more complex than in the 
case of the usual statistical mechanical models, requiring great 
computational effort and specific numerical techniques for the 
production of the data. Moreover, we must consider three lim- 
its in order to obtain the desired physical results from the sim- 
ulation data: 1) the infinite-volume limit (or thermodynamic 
limit), 2) the continuum limit (i.e. the value of a must be suf- 
ficiently small when compared to the relevant distance for the 
problem) and 3) the chiral limit (in order to consider physical 
values for the masses of the light quarks). The above limits 
are not independent, since to get to the continuum limit and 
to be able to consider small masses for the quarks one needs 
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a sufficiently large number of lattice points (corresponding to 
a small enough lattice spacing and to a large enough physical 
size of the lattice), which increases considerably the compu- 
tational effort. 

The study of lattice QCD constitutes a so-called Grand 
Challenge computational problem |4]. Indeed, simulations 
of full QCD — i.e. including effects of dynamical fermions 
— for quark masses in the region of physical values are 
still extremely slow and they are in general carried out 
on supercomputers, involving the effort of large collabora- 
tions such as the UKQCD in the United Kingdom and the 
JLQCD in Japan. Also, several research groups have built 
QCD-dedicated computers, using parallel architecture. Ex- 
amples are the Hitachi/CP-PACS machine at the University 
of Tsukuba in Japan |5|, the QCDSP and QCDOC machines at 
Columbia University in the USA |6], and the APE machines 
0] at various research centers in Europe. These computers 
range from about 1 to 10 teraFLOPS. In addition to these 
large projects, many groups base their simulations on clus- 
ters of workstations or personal computers (PC's) |8]. These 
systems do not yet provide the same efficiency in paralleliza- 
tion as the machines with parallel architecture, but their cost 
is much lower. In addition to the computational power, the 
numerical and analytical techniques used in the simulations 
and in the interpretation of the produced data are of great im- 
portance in the field. Significant progress has been achieved 
through the development of more efficient simulation algo- 
rithms, new methods for interpolation and extrapolation of 
the numerical data and a better understanding of the system- 
atic effects to which the simulation may be subject, such as 
finite-volume effects and discretization errors. Progresses in 
the field are reported annually at the Lattice conference |9]. 

Despite the great computational difficulty, numerical stud- 
ies of QCD have already provided important contributions to 
the study of the strong force 1 10] . In particular, numerical 
simulations of QCD are now able to produce calculations of 
the strong coupling constant u s (po), taken at a fixed reference 
scale /jo, with precision comparable to the experimental one or 
better ll ill . These results are presently included in the world 
average for this quantity 1 12]. Also, the mass spectrum of the 
light hadrons (including the two light quarks and the strange 
quark) has been determined (with great precision) 1 13] for the 
quenched case, in which the configurations are produced with- 
out considering effects of dynamical quarks. One does not 
obtain complete agreement with the experimental spectrum, 
but the observed discrepancies are of at most 10%. Similar 
calculations are now being performed for the full-QCD case. 
Finally, lattice simulations constitute the only known evidence 
for the quark-deconfining transition at finite temperature lll4Tl 
and its predictions are of direct interest for the current exper- 
iments in search of new states of matter in the laboratories 
Brookhaven and CERN. 

There is presently great interest in the results of the simula- 
tions described above and one hopes to be able to solve many 
theoretical questions about QCD and the standard model 
lH5ll . Indeed, lattice-QCD simulations are now able to pro- 
vide quantitative predictions with errors of a few percent. 
This means that these simulations will soon become the main 



source of theoretical results for comparison with experiments 
in high-energy physics 11611 . enabling a much more complete 
understanding of the physics of the strong force. 



A. Lattice QCD at the IFSC-USP 

Since the beginning of 2001 we have been carrying out a 
project on numerical simulations of lattice gauge theories at 
the Physics Department of the University of Sao Paulo in Sao 
Carlos (IFSC-USP), funded by FAPESP 0. The project 
included the installation of 2 PC clusters (with a total of 28 
processing nodes). The resulting computer power is of ap- 
proximately 40 gigaFLOPS for peak performance. We have 
performed production runs since July of 2001 and have started 
intensive parallel simulations in November 2002. Our main 
research topic (see Sections lTIMVlll below) is the investigation 
of the infrared behavior of various propagators and vertices in 
Landau gauge with the goal of verifying the so-called Gribov- 
Zwanziger confinement scenario I18lfl9il . In order to reduce 
the computational cost of the simulations, we consider the 
pure SU(2) gauge theory, including studies in three (instead 
of the usual four) space-time dimensions. 

Besides the topics described below, we also carry out nu- 
merical studies of: gauge-fixing algorithms 12011 . Gribov-copy 
effects 1 21], the chiral phase transition of QCD with two dy- 
namical fermions 1 22], the equation of state of spin models 
with Goldstone modes 12311 - cluster percolation 12411 and short- 
time dynamics for spin models 12511 . 



n. CONFINEMENT SCENARIOS IN LANDAU GAUGE 

As said above, the study of the infra-red (IR) limit of QCD 
is of central importance for understanding the mechanism of 
confinement and the dynamics of partons at low energy. De- 
spite being non-gauge-invariant, gluon and ghost propagators 
are powerful tools in the (non-perturbative) investigation of 
this limit 1 26, 27]. In fact, according to the Gribov-Zwanziger 
and to the Kugo-Ojima [30] confinement scenar- 
ios in Landau gauge, the ghost propagator must show a diver- 
gent behavior in the IR limit — stronger than p~ 2 — for van- 
ishing momentum p. This strong IR divergence corresponds 
to a long-range interaction in real space, which may be re- 
lated to quark confinement. At the same time, according to 
the former scenario, the gluon propagator must be suppressed 
and may go to zero in the IR limit [18, 28ll29il3lll . This would 
imply that the real-space gluon propagator is positive and neg- 
ative in equal measure, i.e. reflection positivity is maximally 
violated 1321 133I1 . As a consequence, the Euclidean 2-point 
function cannot represent the correlations of a physical par- 
ticle. This result may be viewed as an indication of gluon 
confinement |2g|. 

These theoretical predictions have been confirmed by stud- 
ies using Dyson-Schwinger equations (DSE's) J26l E7ll . In 
part icular, studies of DSE's in Landau gauge have found 
10 l2l an IR behavior of the form G(p) ~ p- 2 *- 2 = 
p-2a -2 f or tne gjjogj- propagator and of the form D{p) ~ 
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p-"- - = p<-"D L for the gluon propagator with the same expo 



4k-2 =/ ,2« d -2 ^ 

nent K (i.e. with = ap/l). In 4d one usually finds K > 0.5 
for pure SU(N C ) gauge theory. Note that K > 0.5 implies 
Z)(0) = 0. For the 3d case the exponents are a n ~ 0.4 and 
cid ~ 1.3. Note that in the c/ dimensional case l34l l35ll the 
relation between au and an is given by = 2«g + (4 — cQ/2, 
implying for the quantity 0C s (p) = (g 2 /An)D(p)G 2 (p)p 6 the 
IR behavior p 2 ^ a o- 2a a) — p 4 - rf . Thus, in the 4-dimensional 
case the running coupling OL s (p) displays and IR fixed point. 

Numerical studies of lattice gauge theories confirm the IR 
divergence of the Landau ghost propagator 82 ll 1371 13811 and 
an IR suppression of the gluon propagator. More precisely, a 
decreasing gluon propagator at small momenta has been ob- 
tained for the 3d 51/(2) Landau case using very large lattices 
13^1 and — recently — in the Ad SU{3) Landau case with the 
use of asymmetric lattices |40]. Similar results has also been 
obtained for the the equal-time three-dimensional transverse 
gluon propagator in Ad SU (2) Coulomb gauge |41]. In this 
last case, one also obtains an excellent fit of the transverse 
propagator by a Gribov-like formula. Finally, direct support 
to the Gribov-Zwanziger and to the Kugo-Ojima scenarios has 
been presented in 1 42] and in |43], respectively. 

Thus, the two nonperturbative approaches above seem to 
support the Gribov-Zwanziger and the Kugo-Ojima confine- 
ment scenarios in Landau gauge. However, the agreement be- 
tween the two methods is still at the qualitative level. More- 
over, recent lattice studies l44ll45Tl seem to indicate a null IR 
limit for OL s (p), instead of a finite nonzero value. At the same 
time, a study based on DSE's |46] showed that torus and con- 
tinuum solutions are qualitatively different. This suggests a 
nontrivial relation between studies on compact and on non- 
compact manifolds and could have important implications for 
lattice studies. Also, the nonrenormalizability of the ghost- 
gluon vertex — proven at the perturbative level l47tl . con- 
firmed on the lattice E EH (for p > 200 MeV) and used 
in DSE studies to simplify the coupled set of equations — 
has been recently criticized in Ref. 1 44] . Thus, clear quantita- 
tive understanding of the two confinement scenarios is still an 
open problem. 



III. INFINITE-VOLUME LIMIT 

The study of the IR behavior of propagators and vertices, 
i.e. for momenta smaller than 1 GeV, requires careful consid- 
eration of the infinite-volume limit. Indeed, since the small- 
est non-zero momentum that can be considered on a lattice 
is given by p m i„ ri2k/L — where L is the size of the lattice 
in physical units — it is clear that one needs to simulate at 
very large lattice sizes in order to probe the small-momentum 
limit. The consideration of very large lattice sizes requires 
parallelization and high efficiency of the code in order to ob- 
tain good statistics in the Monte Carlo simulation. Thus, an 
optimized parallel code is of great importance E^l . Our nu- 
merical code is parallelized using MPI; for the random num- 
ber generator we use a double-precision implementation of 
RANLUX (version 2. 1) with luxury level set to 2. 
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FIG. 1: Plot of the rescaled gluon propagator aD(0) /Z(a) as a func- 
tion of the momentum for V = 80 3 and (3 = 4.2 ( x ), 5.0 (□), 6.0 (0). 
Error bars are obtained from propagation of errors. 



A. Very large lattice side 

In Ref. 13^1 we have evaluated the lattice gluon propagator 
D(k) and study it as a function of the momentum p(k) in the 
3d St/ (2) Landau case, using data from the largest lattice side 
to date, i.e. up to 140 3 . This allowed us to consider momenta 
as small as 5 1 MeV (in the deep IR region) and physical lattice 
sides almost as large as 25 fm. 

In order to compare lattice data at different P's, we apply the 
matching technique described in 15011 . i.e. the propagators are 
multiplied by a factor Z(a) depending on p or, equivalently, 
on the lattice spacing a. The method works very well (see 
Fig-Q; indeed, data obtained using different P values nicely 
collapse into a single curve. We find that the gluon propaga- 
tor decreases in the IR limit for momenta smaller than p max , 
which corresponds to the mass scale M in a Gribov-like prop- 
agator D{p) = p 2 1 (p 4 + M 4 ) . From the plot we can estimate 
p max =M = 0.8+j^ Vo = 350+^° MeV, in agreement with 
Ref. jsHl . (Here G is the string tension.) 

In Fig.|2]we plot the rescaled gluon propagator at zero mo- 
mentum, namely aD(Q)/Z(a), as a function of the inverse lat- 
tice side L~ l = 1/ (aN) in physical units (fm -1 ). We see that 
aD(0) /Z(a) decreases monotonically as L increases, in agree- 
ment with Ref. l52ll . It is interesting to notice that these data 
can be well fitted using the simple Ansatz d + b/L c both with 
d = and d ^ (see Figure |2ji. In order to decide for one 
or the other result one should go to significantly larger lattice 
sizes. 

Also, in Ref. l39ll we have shown that the data for the gluon 
propagator are well fitted by Gribov-like formulae, yielding 
an IR critical exponent K « 0.65 in agreement with recent an- 
alytic results (see SectionlllV Recently 1 53] we have extended 
this analysis to the ghost propagator, considering lattice vol- 
umes up to 80 3 for the coupling p = 4.2. A fit to the data 
using the fitting function bj p" (in the interval p < 0.5 GeV) 
gives a = 2.40(2). This result would imply ao ~ 0.4, also in 
agreement with the results reported in Sectionllllabove. 
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B. Asymmetric lattices 

Recently, very asymmetric lattices 1I40I l54ll have been con- 
sidered in order to explore the IR limit od QCD. As a test 
of this method, we have extended the gluon propagator 
study presented in 13911 [for the 3d ST/ (2) case in minimal 
Landau gauge], by including results for the ghost propagators 
from very large lattices. At the same time, we evaluated the 
propagators using also asymmetric lattices, in order to verify 
possible systematic effects related to the use of asymmetric 
lattices (as suggested in 14510 . by comparing the results to the 
ones obtained for symmetric lattices. 

We find, for both propagators, clear evidences of system- 
atic effects at relatively small momenta, i.e. p < 1.5^0 rs 
650MeV. In particular, the gluon (respectively, ghost) prop- 
agator is less suppressed (respec. enhanced) in the IR limit 
when considering asymmetric lattices than for the case of 
symmetric lattices (see Fig. [3}. This implies that the esti- 
mates for the IR critical exponents oq and ad are systemat- 
ically smaller in the asymmetric case compared to the sym- 
metric one. 

Also, for the gluon propagator and considering the asym- 
metric lattices, one would estimate a value M < 0.25/a « 0.22 
GeV as a turnover point in the IR, i.e. the momentum p max = 
M for which the propagator reaches its peak. On the other 
hand, considering the largest (symmetric) lattice volumes, i.e. 
V = 140 3 (see the top plot in Fig.[3J, the gluon propagator is 
clearly a decreasing function of p for p < 0.5/a, correspond- 
ing to M < 0.435 GeV. This is in agreement with the result 
reported in Section ITlI Al above. We thus see a difference of 
almost a factor 2 between the momentum- turnover point in the 
symmetric and asymmetric cases. 

Finally, we have seen that the extrapolation to infinite vol- 
ume of results obtained using asymmetric lattices is also most 
likely affected by systematic errors. We conclude that, even 
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FIG. 3: Plot of the gluon propagator D(p) (top) and of the ghost 
propagator G(p) (bottom) as a function of p for lattice volumes V = 
8 2 x 140 (x), 12 2 x 140(D), 16 2 x 140(0) and V = 140 3 (O)- All 
quantities are in lattice units. Note the logarithmic scale on the y axis 
in the bottom plot. Errors represent one standard deviation. 




1 ' 1 1 1 1 

0.05 0.1 0.15 0.2 

1/L (fm" 1 ) 

FIG. 2: Plot of the rescaled gluon propagator at zero momentum as 
a function of the inverse lattice side. We also show the fit of the data 
using the Ansatz d + b/L c , both with d = and with d ^= 0. Error 
bars are obtained from propagation of errors. 



though using an asymmetric lattice does not modify the qual- 
itative behavior of the two propagators, one should be careful 
in extracting quantitative information from such studies. 

IV. REFLECTION-POSITIVITY VIOLATION 

The relation between reflection positivity and Euclidean 
correlation functions can be made explicit by considering the 
spectral representation 1 2fiJ B^l 

D(p) = rdn?4p- 2 (1) 
Jo p M -\-mr 

for the Euclidean propagator in momentum space. Then, the 
statement of reflection positivity is equivalent to a positive 
spectral density p(f« 2 ). This implies that the temporal cor- 
relator at zero spatial momentum D(t ,p = 0) can be written 

as 

C(t) = D(t,0) = f d(op{(0 2 )e- at . (2) 
Jo 
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FIG. 4: Real-space propagator C(f) as a function of / for coupling 
P = 5.0 and lattice volume V = 140 3 . Errors have been evaluated 
using bootstrap with 1000 samples. All quantities are in lattice units. 



FIG. 5: Fit of C(t) as a function of t (in fm) using a sum of two 
Gribov-like propagators for lattice volume V = 140 3 and coupling 
(3 = 4.2. We also display the two Gribov-like propagators separately. 



Clearly, a positive density p(tt> 2 ) implies that C(f) > 0. No- 
tice that having C(t) positive for all t does not ensure the posi- 
tivity of p((0 2 ). On the other hand, finding C(t) < for some t 
implies that p(tt> 2 ) cannot be positive, suggesting confinement 
for the corresponding particle. 

On the lattice, the real-space propagator can be evaluated 
using 
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-2k ik t/N 



D(k ,0) 



(3) 



where N is the number of points per lattice side and D(k) is 
the propagator in momentum space. If the lattice action sat- 
isfies reflection positivity |57], then we can write the spectral 
representation 

-E„t 



(4) 



438MeV = 0.995\/o, where a is the string tension. Let us 
also recall that the momentum-space Gribov-like propaga- 
tor has its maximum at p max = M (see Section ITlI Al above). 
Moreover, we find that finite-size effects seem to become 
important only at t > 3 /to. This means that our data for 
t £ [0,3] fm are essentially infinite-volume continuum results. 
In the scaling region, the data are well described by a sum 
of Gribov-like formulas (see Fig. [5}, with a light-mass scale 
tow 0.74^0 = 325 MeV. 

It has been suggested 1 60, 61] that the violation of spectral 
positivity in lattice Landau gauge be related to the quenched 
auxiliary fields used for gauge fixing. We note that the fitting 
form proposed for C(t) in [61] describes reasonably well our 
data up to t = 3 fm — yielding a light-mass scale of about 
1.14^/a = 500 MeV — but cannot account for the oscillatory 
behavior observed at very large separations. 



where r„ are positive-definite constants. Clearly, this implies 
that C(t) is non-negative for all values of t . 

Numerical indications of a negative real-space lattice Lan- 
dau gluon propagator have been presented in the 3d 51/(2) 
case 15811 . in the magnetic sector of the 4-d SI/ (2) case at finite 
temperature [59] and, recently, in the Ad SI/ (3) case for one 
"exceptional" configuration |33[. 

Using data from the largest lattice sides to date, we ver- 
ify explicitly (in the 3d case) |32] the violation of reflection 
positivity for the SI/ (2) lattice Landau gluon propagator (see 
Fig. |4}. In particular, the propagator becomes negative at 
t ps 0.7 fm and the minimum is reached at f m ,-„ pa 1/to (see 
Fig.|5}- Note that the Gribov-like propagator 
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has its minimum at f„„„ = 7t/(MV2). Thus, the above re- 
sult for t m in would imply M pa %j \flfm~ l pa 2.22 /m _1 = 



V. GHOST-GLUON VERTEX 

In the framework of quantum field theory, Faddeev-Popov 
ghosts are introduced in order to quantize non-Abelian gauge 
theories. Although the ghosts are a mathematical artifact and 
are absent from the physical spectrum, one can use the ghost- 
gluon vertex and the ghost propagator to calculate physical 
observables, such as the QCD running coupling OL s (p), using 
the relation 



a s (p) = oto 



Z3(p)Z 2 3 (p) 



(6) 



Here oto = go/An is the bare coupling constant and Z^{p), 
Zi(p) and Z\ (p) are, respectively, the gluon, ghost and ghost- 
gluon vertex renormalization functions. The above formula 
gets simplified if one considers the Landau gauge. Indeed, 
in this case the vertex renormalization function Z\ (p) is finite 
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FIG. 7: The scalar quantity G [q,k,§), defined in Eq. (27) of Ref. 
ItjBll . as a function of the magnitude of the incoming anti-ghost mo- 
mentum q. Full symbols correspond to P = 4.2 and open symbols to 
FIG. 6: Results for Z\ (p) for the lattice volume V = 16 4 as a p = 6.0; circles are used for V = 20 3 and triangles for V = 30 3 . We 



function of p = p/a in GeV, considering symmetric momenta, i.e. 
with 4 equal components. We show data obtained using two differ- 
ent gauge-fixing methods (with and without the so-called smearing 
method 1 64]). Error bars were evaluated using the bootstrap method 
with 250 samples. 



and constant, i.e. independent of the renormalization scale p, 
at least to all orders of perturbation theory iEtH^ . Of course, 
it is important to verify in a non-perturbative way that this re- 
sult really holds. If it does, one can consider in the Landau 
gauge a definition of the running coupling constant that re- 
quires only the calculation of the gluon and ghost propagators 

In Refs. j48il we have studied the reduced ghost-gluon ver- 
tex function T? bc (0,p) and the renormalization function Z\ (p) 
in minimal Landau gauge at the asymmetric point (Q;p,—p) 
in the St/ (2) case. We find that the vertex function has the 
same momentum dependence of the (lattice) tree-level vertex 
— i.e. ~ p cos (71 pa/L) ~ sm(2%pa/ L) — and that Z\(p) is 
approximately constant and equal to 1, at least for momenta 
p > 1 GeV (see Fig.|6}. This is a direct non-perturbative veri- 
fication of the perturbative result reported above. In particular, 
using the result obtained at the largest value of p considered 
(i.e. P = 2.4) we obtained Z^ l (p) = 1.02^, where errors in- 
clude Gribov-copy effects and discretization errors related to 
the breaking of rotational invariance. 

Recently, this study has been extended (in the 3d case) |65] 
to other kinematical configurations, including in particular the 
symmetric point p 2 = q 2 = k 2 . We find that the vertex is es- 
sentially constant and of order one, for all momentum con- 
figurations (see Fig.Q, confirming the results obtained in the 
4d case 1451 14811 . This result is also in agreement with pre- 
dictions from functional methods in the 3d case |66]. In the 
same reference we also present the first numerical study of the 
Landau-gauge three-gluon vertex in the 3d case and results for 
the spectrum of the Faddeev-Popov operator, which plays an 
important role in the Gribov-Zwanziger scenario fT8ll28Ll67ll . 



plot data for the case with the three momenta equal. 



In particular we have shown that (see Fig. 7 in Ref. |65]) 
the smallest non-zero eigenvalue of the Faddeev-Popov ma- 
trix goes to zero when the infinite-volume limit is approached. 
As a consequence, in the continuum limit, the average lattice 
Landau configuration should belong to the first Gribov hori- 
zon, supporting the Gribov-Zwanziger mechanism of confine- 
ment (see also 1 42]). 



VI. GHOST CONDENSATES 

The QCD vacuum is known to be highly non-trivial at 
low energies 16811 . This non-trivial structure manifests itself 
through the appearance of vacuum condensates, i.e. vacuum 
expectation values of certain local operators. In perturba- 
tion theory these condensates vanish, but in the SVZ-sum- 
rule approach 1 69] they are included as a parametrization of 
non-perturbative effects in the evaluation of phenomenologi- 
cal quantities. The two main such operators are O^F^F^ and 
mqyyjtyq, their vacuum expectation values are the so-called 
gluon and quark condensates. Both of these operators have 
mass dimension four. 

In recent years, (gauge-dependent) condensates of mass di- 
mension two have also received considerable attention l29l 
EH S III Eli El El- In par- 
ticular, the gauge condensate (A^A^) has been largely stud- 
ied, since it should be sensitive to topological structures such 
as monopoles 1 70] and it could play an important role in 
the quark-confinement scenario through monopole conden- 
sation 1 83]. Moreover, the existence of a gauge condensate 
would imply a dynamical mass generation for the gluon and 
ghost fields 129. 17TI1 . Possible effects of the gauge conden- 
sate (A*A*) on propagators and vertices (in Landau gauge) 
have been studied through lattice simulations in Refs. 1 44, 84], 
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yielding (a£a£) « 3 GeV 2 . 

Other vacuum condensates of mass dimension two con- 
sidered by several groups are the ghost condensates. These 
condensates were first introduced in SU(2) gauge t heory in 
maximally Abelian gauge (MAG) S 0,^3]. More 
recently, the same condensates have been studied in other 
gauges ins, such as the Curci-Ferrari and the Landau 
gauges. In all cases it was found that the ghost condensates 
are related to the breakdown of a global SL(2,R) symmetry 
j2^,E3- 1° MAG the diagonal and off-diagonal components 
of the ghost propagators are modified 1 72[ 17711 by ghost con- 
densation. Similar results have been obtained in other gauges 
(z3,[z3]- I n particular, in Landau gauge it was shown |78] that 
the off-diagonal (anti-symmetric) components of the ghost 
propagator G cd {p) are proportional to the ghost condensate 
v. 

Finally, a mixed gluon-ghost condensate of mass dimension 
two has also been studied by several authors ll73l l75l l79l l80l 
lilll . using various gauges. This mixed condensate is of partic- 
ular interest when considering interpolating gauges |8(J. In- 
deed, it allows to generalize and relate results obtained in dif- 
ferent gauges for the gauge condensate (AjjA*) and the ghost 
condensates. Moreover, in MAG this mixed condensate would 
induce a dynamic mass for the off-diagonal gluons |8lJ, giv- 
ing support to the Abelian-dominance scenario 1 86]. Thus, the 
various gauge and ghost condensates could all play an impor- 
tant role in the dual superconducting scenario of quark con- 
finement l87ll . being related to monopole condensation and to 
Abelian dominance. 

In Ref. 18811 we carried out a thorough investigation of ghost 




FIG. 8: Results for the Binder cumulant U for the quantity ty b (p) as 
a function of the lattice side L (in GeV~ 1 ) for various lattice volumes 
V and momentum with p=N/4, We show the data corresponding to 
asymmetric momenta [for (3 = 2.2 (□), (3 = 2.3 (O) and P = 2.4 (A)] 
and to symmetric momenta (with the corresponding filled symbols 
for each (3). Errors have been estimated using the bootstrap method 
with 10,000 samples. 



condensation in the so-called Overhauser channel for pure 
ST/ (2) Yang-Mills theory in minimal Landau gauge. In par- 
ticular, we evaluate numerically the off-diagonal components 
of the ghost propagator G cd (p) as a function of the momen- 
tum p. We find that (§ b {p)) = E bcd (G cd (p)) /2 is zero within 
error bars, but with large fluctuations. At the same time, we 
see clear signs of spontaneous breaking of a global symmetry, 
using the quantity § h (p) as an order parameter. As in the case 
of continuous-spin models in the ordered phase (see for exam- 
ple 18*^0 . spontaneous symmetry breaking is supported by two 
(related) observations: 1) by comparing the statistical fluctu- 
ations for the quantities § b (p) and 2) from the non- 
Gaussian shape of the statistical distribution of § b {p), which 
can be observed by considering a histogram of the data or by 
evaluating the so-called Binder cumulant (see Fig. [8}. Since, 
in Landau gauge, the vacuum expectation value of the quan- 
tity § b (p) should be proportional (in the Overhauser channel) 
to the ghost condensate v 178], it seems reasonable to conclude 
that the broken symmetry is the SL(2,R) symmetry, which is 
related to ghost condensation [ 26, 76] . This interpretation has 
been recently criticized in Ref. 1 90] . There, the nonzero value 
obtained for the Binder cumulant has been explained by con- 
sidering multi-dimensional Gaussian distributions and a mod- 
ified definition for the Binder cumulant. We note however 
that the order parameter considered in our study, i.e. the mag- 
nitude of § b {p), is a scalar quantity. Therefore, in an alog y 
with studies of (9(A^)-vector models (see for example 19 IB ), 
the standard definition of the Binder cumulant should apply to 
our study as well. 

Let us note that, from our data, the Binder cumulant U 
seems to be approximately null at small lattice volume and to 
converge to a value U w 0.45 for physical lattice side L > 15 
Gev -1 w 3 fm, corresponding to a mass scale of less than 100 
MeV. 

In Ref. llssll we have also shown that the sign of § b (p) is 
related to the sign of the Fourier-transformed gluon field com- 
ponents A (q) and that § b (p) has discretization effects similar 
to those obtained for the ghost-gluon vertex |48]. Then, us- 
ing the rescaled quantity \L 2 § h (p)/ cos(npa/L) |, we find (at 
small momenta) a behavior p~ z with z ~ 4, in agreement with 
analytic predictions 1 78 ] . On the other hand, from our fits 
we find that the ghost condensate v is consistent with zero 
within error bars, i.e. the quantity \L 2 § h (p)/cos(npa/L) \ 
does not approach a finite limit at small momenta, at least for 
p > 0.245 GeV. Using the Ansatz 



G cd { P ) 



p 2 8 cd + ve cd 



(7) 



we obtain for the ghost condensate the upper bound v <C 0.058 
GeV 2 . More precisely, our data rule out values of v greater 
than 0.058 GeV 2 w (240 MeV) 2 but would still be consis- 
tent with a ghost condensate v < 0.025 GeV 2 « (160 MeV) 2 . 
Let us note that, in analytic studies l73l l78l I82L 1921 one finds 
that the ghost condensate induces a tachyonic gluon mass pro- 
portional to y/v, which modifies the dynamic gluon mass re- 
lated to the gauge condensate (A^A b ). Thus, one should ex- 
pect a relatively small ghost condensate v in order to obtain 
a global (non-tachyonic) gluon mass. Let us recall that a dy- 
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FIG. 9: Fit for the running coupling a s (p) using eq. J 101 with C = 
0.072(8), a = 1.9(3), A = 1.31(1) and m = 1.0(6). 



namic gluon mass of the order of a few hundred Me V has been 
considered in several phenomenological studies |93]. A sim- 
ilar mass scale was also obtained in numerical studies of the 
gluon propagator in Landau gauge ll32ll3^.l5lll . 

VII. RUNNING COUPLING 

Of great importance for phenomenological purposes is the 
running coupling strength a s (p) 16311 defined in Eq. ijfji above. 
In particular, this quantity enters directly the quark Dyson- 
Schwinger equation (DSE) and can be interpreted as an effec- 
tive interaction strength between quarks 1 94] . Let us note that, 
working in Landau gauge and in the momentum-subtraction 
scheme, the running coupling can be written as |63] 

a s ( P ) = aoF(p)J 2 (p) , (8) 

where F(p) and J(p) are, respectively, the gluon and the ghost 
form factors and we used the result Z\ (p) = 1. As explained 
in Sec.|II|above, studies using DSE's have found that, if the IR 
sum rule 2ciq — ao = is satisfied, then this running coupling 
develops a fixed point in the IR limit (see for example 1 35]) 

limoCs(p) = OCc = constant, (9) 

p->0 

with a c « S.92/N c in the SU(iV c ) case for K w 0.596. 

This quantity has been studied numerically by several 
groups (see for example 1331 l44l l95lo . In Ref. ll96ll we have 
evaluated this running coupling constant and tried a fit to the 
data using the fitting function 

a s (p) =Cp 4 /[(p 4 +m)s(a)] (10) 

where s(a) = (ll/247t 2 )log [1 + (p 2 /A 2 )"}. Note that if a = 
2 one finds a fixed point in the IR limit equal to a c = 



P 

FIG. 10: Plot of the gluon propagator D(p) as a function of p for 
lattice volume V = 52 4 at (3 = 2.2. All quantities are in lattice units 
(here a ~ 1 GeV -1 ). Errors represent one standard deviation. The 
larger value at p = is probably due to an insufficient number of 
gauge-fixing iterations. 

247t 2 CA 4 /l lm. Also, this fitting functions satisfies the lead- 
ing ultraviolet behavior ~ 1/log [(p 2 /A 2 )] for the running 
coupling. Results are reported in Fig.|9j we find a = 1.9(3) 
and 0C C = 24tt. 2 CA 2 /1 lm w 4.6, in agreement with the result 
reported above. Similar results have also been obtained in 
H38I1 . We stress that in Refs. 1381 19611 we have compared re- 
sults obtained using two slightly different lattice formulations, 
yielding consistent results in all cases considered. 

As said in Section UTI above, recent lattice studies J44ll45ll 
seem to indicate a null IR limit for OL s (p), instead of a fi- 
nite nonzero value a c . We should stress, however, that in 
these studies special care must be taken in order to eliminate 
finite-size effects, especially when the IR region is consid- 
ered. Indeed, we know that, if the the gluon propagator is 
not suppressed at small momenta, the true infrared regime is 
surely not reached yet. This is a difficult numerical task since 
we have checked that in the 4d ST/ (2) case a lattice volume 
V = 52 4 at p = 2.2, corresponding to a physical lattice side of 
about 10 fm, is still not sufficient to show a decreasing prop- 
agator in the limit of small momenta (see Fig. I10i . We note 
that in Fig. 1 of [44] and in Fig. 3 of |45] the running cou- 
pling starts to decrease for a momentum p « 0.5 GeV. Thus, 
from our discussion in Section El Al above. considering the re- 
sults shown in Fig.[H)]and also the results reported in 
we can say that the data considered in |44, 45] for momenta 
p < 500 MeV are most likely affected by strong finite-size 
effects. 



VIII. CONCLUSIONS 

The Gribov-Zwanziger and the Kugo-Ojima confinement 
scenarios in Landau gauge are well supported, at the qualita- 
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tive level, by several studies based on different analytical and 
numerical approaches. On the other hand, as we have shown 
here, a comprehensive analysis of these confinement scenar- 
ios at the quantitative level using lattice numerical simulations 
could represent a very challenging task. 



tado de Sao Paulo (FAPESP) (under grant # 00/05047-5). 
Partial support from Conselho Nacional de Desenvolvimento 
Cientffico e Tecnologico (CNPq) is also acknowledged. 



ACKNOWLEDGMENTS 



The authors thank the organizers for the invitation. Re- 
search supported by Fundacao de Amparo a Pesquisa do Es- 



[1] An Elementary Primer for Gauge Theory, K. Moriyasu (World 
Scientific, Singapore, 1983). 

[2] K.G. Wilson, Phys. Rev. D 10, 2445 (1974); Lattice gauge the- 
ories. An introduction, H.J. Rothe (World Scientific, Singapore, 
1997). 

[3] A guide to Monte Carlo simulations in Statistical Physics, D.P. 

Landau and K. Binder (Cambridge University Press, Cam- 
bridge UK, 2000). 
[4] Parallel computer architecture: a hardware/software approach, 

D. E. Culler, J.P. Singh, and A. Gupta (Morgan Kaufmann, San 

Francisco CA, 1999). 
[5] K. Kanaya, Fortsch. Phys. 50, 531 (2002). 
[6] R.D. Mawhinney, Parallel Comput. 25, 1281 (1999); PA. Boyle 

et al, J. Phys. Conf. Sen 16, 129 (2005). 
[7] A. Bartoloni et al, Nucl. Phys. Proc. Suppl. A 60, 237 (1998); 

R. Ammendola et al, Nucl. Phys. Proc. Suppl. B 119, 1038 

(2003). 

[8] M. Liischer, Nucl. Phys. Proc. Suppl. B 106, 21 (2002). 
[9] Proceedings of the XXIIIrd International Symposium on Lattice 
Field Theory, PoS LAT2005 and previous years. 
[10] M. Creutz, AIP Conf. Proc. 690, 52 (2003). 
[11] A. Cucchieri, in the Proceedings of Hadron Physics 2002, 
edited by C.A.Z. Vasconcellos et al. (World Scientific, Cinga- 
pura, 2003) 137. 

[12] Quantum chromodynamics, I. Hinchliffe, in Review of particle 
physics, Particle Data Group (K. Hagiwara et al), Phys. Rev. D 
66(2002)010001. 

[13] S. Aoki etal, Phys. Rev. Lett. 84, 238 (2000). 

[14] F. Karsch, Lect. Notes Phys. 583, 209 (2002). 

[15] F. Wilczek, Nucl. Phys. Proc. Suppl. B 119, 3 (2003). 

[16] C.T.H. Davies et al, Phys. Rev. Lett. 92, 022001 (2004); P. 
Lepage and C. Davies, Int. J. Mod. Phys. A 19, 877 (2004). 

[17] |http : / /lattice ■ if sc .usp .br/| 

[18] V.N. Gribov. Nucl. Phys. B 139, 1 (1978). 

[19] D. Zwanziger. Nucl. Phys. B 412, 657 (1994). 

[20] A. Cucchieri and T. Mendes, Nucl. Phys. B 471, 263 (1996); 
Nucl. Phys. Proc. Suppl. B 53, 811 (1997); Phys. Rev. D 57, 
3822(R) (1998); Comput. Phys. Commun. 154, 1 (2003). 

[21] A. Cucchieri, Nucl. Phys. B 508, 353 (1997). 

[22] T. Mendes and A. Cucchieri, AIP Conf. Proc. 739, 599 (2004); 
T. Mendes, Nucl. Phys. A 702, P29 (2002); in the Proceedings 
of Hadron Physics 2002, edited by C.A.Z. Vasconcellos et al. 
(World Scientific, Cingapura, 2003) 297; AIP Conf. Proc. 756, 
413 (2005). 

[23] A. Cucchieri and T. Mendes, J. Phys. A 38, 4561 (2005); Nu- 
merical simulation of 'N -vector spin models in a magnetic field, 
to appear in Braz. J. Phys. 



[24] W.G. Wanzeller, A. Cucchieri, G. Krein and T. Mendes, Braz. 

J. Phys. 34, 247 (2004). 
[25] Dynamic behavior of cluster observables for the 2d Ising 

model, W.G. Wanzeller, T. Mendes and G. Krein, to appear 

in Braz. J. Phys; W.G. Wanzeller, T. Mendes and G. Krein, 

Icond-mat 70 60543*51 
[26] R. Alkofer and L. von Smekal, Phys. Rept. 353, 281 (2001). 
[27] A. Holl, CD. Roberts and S.V. Wright, nucl-th/ 0601071. 
[28] D. Zwanziger, Nucl. Phys. B 412, 657 (1994). 
[29] R.F Sobreiro, S.P Sorella, D. Dudal and H. Verschelde, Phys. 

Lett. B 590, 265 (2004). 
[30] T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66, 1 (1979); 

Erratum - Prog. Theor. Phys. 71, 1121 (1984); T. Kugo, 

|hep-th/9511033| 
[31] D. Zwanziger, Phys. Lett. B 257, 168 (1991); Nucl. Phys. B 

364, 127 (1991); Nucl. Phys. B 378, 525 (1992). 
[32] A. Cucchieri, T. Mendes and A.R. Taurines, Phys. Rev. D 71, 

051902(R) (2005) 
[33] S. Furui and H. Nakajima, Phys. Rev. D 70, 094504 (2004). 
[34] D. Zwanziger, Phys. Rev. D 65, 094039 (2002). 
[35] C. Lerche and L. von Smekal, Phys. Rev. D 65, 125006 (2002). 
[36] R. Alkofer, W. Detmold, C.S. Fischer and P. Maris, Phys. Rev. 

D 70, 014014 (2004). 
[37] H. Suman and K. Schilling, Phys. Lett. B 373, 314 (1996); S. 

Furui and H. Nakajima, Phys. Rev. D 69, 074505 (2004). 
[38] J.C.R. Bloch, A. Cucchieri, K. Langfeld and T. Mendes, Nucl. 

Phys. B 687, 76 (2004). 
[39] A. Cucchieri, T. Mendes and A.R. Taurines, Phys. Rev. D 67, 

091502(R) (2003). 
[40] O. Oliveira and P.J. Silva, AIP Conf. Proc. 756, 290 (2005); PoS 

LAT2005, 286 (2005); PoS LAT2005, 287 (2005); hep-lat 

/0511043. 

[41] A. Cucchieri and D. Zwanziger, Phys. Rev. D 65, 014001 

(2002); Phys. Lett. B 524, 123 (2002). 
[42] A. Cucchieri, Nucl. Phys. B 521, 365 (1998). 
[43] H. Nakajima, S. Furui and A. Yamaguchi, hep-lat/0007001 

S. Furui and H. Nakaiima, |hep-lat/040302l| Phys. Rev. D73, 

074503 (2006). 

[44] P. Boucaud et a/., |hep-ph/0507104| 

[45] E.M. Ilgenfritz, M. Muller-Preussker, A. Sternbeck and A. 

Schiller, |hep-lat/0601027| 
[46] C.S. Fischer, B. Gruter and R. Alkofer, hep-ph/0506053 
[47] J.C. Taylor, Nucl. Phys. B 33, 436 (1971). 
[48] A. Cucchieri, T. Mendes and A. Mihara, JHEP 0412, 012 

(2004); AIP Conf. Proc. 739, 602 (2005). 
[49] A. Cucchieri, T. Mendes, G. Travieso and A.R. Taurines, 

|hep-lat 7030 BOO? in the Proceedings of the 15th Symposium 



9 



on Computer Architecture and High Performance Computing, 

edited by L.M. Sato et al. (IEEE Computer Society Press, Los 

AlamitosCA, 2003) 123. 
[50] D.B. Leinweber et al, Phys. Rev. D 60, 094507 (1999) [Erra- 

tum-ibid. D 61, 079901 (2000)]. 
[51] A. Cucchieri, Phys. Rev. D 60, 034508 (1999). 
[52] ED. Bonnet et al., Phys. Rev. D 64, 034501 (2001). 
[53] A. Cucchieri, A. Maas and T. Mendes, in preparation. 
[54] M.B. Parappilly et a/., ihep-lat/060ToTol 
[55] A. Cucchieri and T. Mendes, Phys. Rev. D 73, 071502(R) 

(2006). 

[56] H. Aiso et al, Nucl. Phys. Proc. Suppl. B 53, 570 (1997). 

[57] Quantum Fields on a Lattice, I. Montvay and G. Miinster (Cam- 
bridge University Press, Cambridge UK, 1994). 

[58] A. Cucchieri, F. Karsch and P. Petreczky, Phys. Lett. B 497, 80 
(2001); Nucl. Phys. Proc. Suppl. B 94, 385 (2001). 

[59] A. Cucchieri, F. Karsch and P. Petreczky, Phys. Rev. D 64, 
036001 (2001). 

[60] C.A. Aubin and M.C. Ogilvie, Phys. Lett. B 570, 59 (2003); 
[61] C. Aubin and M.C. Ogilvie, Phys. Rev. D 70, 074514 (2004). 
[62] O. Piguet and S.P. Sorella, Lect. Notes Phys. M 28 (1995) 1. 
[63] L. von Smekal, R. Alkofer and A. Hauck, Phys. Rev. Lett. 79, 

3591 (1997); Annals Phys. 267 (1998) 1 [Erratum-ibid. 269 

(1998) 182]. 

[64] J.E. Hetrick and P. de Forcrand, Nucl. Phys. Proc. Suppl. B 63, 

838 (1998). 

[65] A. Cucchieri, A. Maas and T. Mendes, |hep-lat/0605011| 
[66] W. Schleifenbaum, A. Maas, J. Wambach and R. Alkofer, Phys. 

Rev. D 72, 014017 (2005); W. Schleifenbaum, M. Leder and H. 

Reinhardt, |hep-th/0605115| 
[67] A. Sternbeck, E. M. Ilgenfritz and M. Muller-Preussker, Phys. 

Rev. D 73, 014502 (2006). 
[68] The QCD vacuum, hadrons and the superdense matter, E.V. 

Shuryak (World Scientific, Singapore, 1988). 
[69] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. 

B 147, 385 (1979); Nucl. Phys. B 147, 448 (1979). 
[70] FV. Gubarev, L. Stodolsky and V.I. Zakharov, Phys. Rev. Lett. 

86, 2220 (2001). 

[71] K.I. Kondo, Phys. Lett. B 560, 44 (2003); D. Dudal etal, JHEP 
0401, 044 (2004); D. Dudal, R.F Sobreiro, S.P. Sorella and H. 
Verschelde, Phys. Rev. D 72, 014016 (2005). 

[72] K.I. Kondo and T. Shinohara, Phys. Lett. B 491, 263 (2000). 

[73] D. Dudal and H. Verschelde, J. Phys. A 36, 8507 (2003). 

[74] K.I. Kondo, Phys. Lett. B 572, 210 (2003). 

[75] K.I. Kondo, T. Murakami, T. Shinohara and T. Imai, Phys. Rev. 
D 65, 085034 (2002). 

[76] D. Dudal et al, JHEP 0212, 008 (2002). 

[77] V.E.R. Lemes, M.S. Sarandy and S.P. Sorella, J. Phys. A 36, 
7211 (2003). 



[78] D. Dudal et al, JHEP 0306, 003 (2003). 

[79] K.I. Kondo, Phys. Lett. B 514, 335 (2001); D. Dudal et al, 
Annals Phys. 308, 62 (2003). 

[80] D. Dudal et al, Phys. Rev. D 70, 114038 (2004). 

[81] D. Dudal et flZ., |hep-th/050l227| 

[82] M.A.L. Capri et al, Phys. Rev. D 73, 014001 (2006). 

[83] G. 't Hooft, Nucl. Phys. B 190, 455 (1981); T. Suzuki, Nucl. 
Phys. Proc. Suppl. B 30, 176 (1993). 

[84] P. Boucaud et al, Phys. Rev. D 63, 114003 (2001); Phys. 
Rev. D 66, 034504 (2002); Phys. Rev. D 67, 074027 (2003); 
|hep-lat/0504017| E. Ruiz Arriola, P.O. Bowman and W. Bro- 
niowski, Phys. Rev. D 70, 097505 (2004); S. Furui and H. 
Nakajima, |hep-lat / 5 3 29] 

[85] M. Schaden HEep^rT7W!??TT]|hep-th/01(18054| in the Pro- 
ceedings of 4th International Conference on Quark confinement 
and the hadron spectrum, edited by W. Lucha and K.M. Maung 
(World Scientific, Singapore, 2001) 258. 

[86] Z.F. Ezawa and A. Iwazaki, Phys. Rev. D 25, 2681 (1982); T. 
Suzuki, Prog. Theor. Phys. 81, 752 (1989). 

[87] Y. Nambu, Phys. Rev. D 10, 4262 (1974); S. Mandelstam, Phys. 
Rept. 23, 245 (1976). 

[88] A. Cucchieri, T. Mendes and A. Mihara, Phys. Rev. D 72, 
094505 (2005). 

[89] Monte Carlo Simulation in Statistical Physics An Introduction, 
K. Binder and D.W. Heermann (Springer- Verlag, Berlin, 2002); 
Computer Studies of Phase Transitions and Critical Phenom- 
ena, O.G. Mouritsen (Springer- Verlag, Berlin, 1984). 

[90] S. Furui and H. Nakajima, |hep-lat / 060"2027l 

[91] K. Kanaya and S. Kaya, Phys. Rev. D 51, 2404 (1995); J. Engels 
and T. Mendes, Nucl. Phys. B 572, 289 (2000). 

[92] H. Sawayanagi, Phys. Rev. D 67, 045002 (2003). 

[93] G. Parisi and R. Petronzio, Phys. Lett. B 94, 51 (1980); F. 
Halzen, G.I. Krein and A.A. Natale, Phys. Rev. D 47, 295 
(1993); M.B. Gay Ducati, F. Halzen and A.A. Natale, Phys. 
Rev. D 48, 2324 (1993); J.H. Field, Phys. Rev. D 66, 013013 
(2002); M.B. Gay Ducati and W.K. Sauter, |r^p^ph/0501057l 

[94] J.C.R. Bloch, Phys. Rev. D 66, 034032 (2002). 

[95] A. Cucchieri, T. Mendes and D. Zwanziger, Nucl. Phys. Proc. 
Suppl. B 106, 697 (2002); K. Langfeld et a/., |hep-th/0209173| 
in the Proceedings of Gargnano 2002, Quark confinement and 
the hadron spectrum, edited by N. Brambilla and G. M. Pros- 
peri (World Scientific, River Edge, 2003) 297; H. Nakajima and 
S. Furui, Nucl. Phys. Proc. Suppl. 129 (2004) 730; A. Stern- 
beck, E.M. Ilgenfritz, M. Muller-Preussker and A. Schiller, PoS 
LAT2005, 333 (2006). 

[96] J.C.R. Bloch, A. Cucchieri, K. Langfeld and T. Mendes, Nucl. 
Phys. Proc. Suppl. B 119, 736 (2003). 



